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Abstract 

We present a systematic theory of Coulomb interaction effects in the nonlinear 
optical processes in semiconductors using a perturbation series in the exciting 
laser field. The third-order dynamical response consists of phase-space filling 
correction, mean-field exciton-exciton interaction, and two-exciton correlation 
effects expressed as a force-force correlation function. The theory provides 
a unified description of effects of bound and unbound biexcitons, including 
memory-effects beyond the Markovian approximation. In the degenerate four- 
wave- mixing experiments, correlation effects are shown leading to polarization 
mixing, ringing, etc. The strong interaction, nonperturbative theory of the 
correlation function is numerically evaluated for a one-dimensional model. 
Approximations for the correlation function are presented. 

71.35.+Z, 42.50.Md 
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I. INTRODUCTION 



Transient four-wave-mixing (FWM) experiments have proven to be a powerful tool in 
probing and understanding optical coherence in semiconductors Subpicosecond spec- 

troscopy yields information on the very early stages of time development of the carrier 
dynamics and many-particle correlations. 

The essential physical picture behind the dynamical evolution of optically excited elec- 
trons and holes can be understood in simple terms []I],|]-|8|. First, the ultrafast dynamics of 
the exciting laser field with frequency near the fundamental band gap of a semiconductor 
creates coherent electron-hole (eh) pairs. Secondly, the motion of the carriers, dominated by 
the Coulomb interaction among them, leads to an ultrafast electric polarization as a source 
of light which can be observed. The scattering of carriers by other carriers, by phonons, and 
by defects leads to polarization decay and loss of optical coherence ||. 

The density-matrix equations of motion which described the dynamics of the micro- 
scopic polarization and particle distribution functions were established by a number of 
groups |T0|-|l4|. Within the mean- field approximation, results for the ultrafast dynamics 
of the electron-hole pairs agreed well with the extant experimental findings, for example, 
for the dynamical Stark effect of the excitons in semiconductors [l5j in a wide range of 
semiconductor bulk and quantum well systems. With recent advances in ultrafast nonlin- 
ear optical spectroscopy and in fabrication of semiconductor heterostructures, the use of 
three-pulse FWM, polarization JTo] , |iT| and phase sensitive measurements of the nonlinear 
polarization [[OJ, and non-degenerate transient FWM ]19|] has led to effects beyond the 
mean-field approximation for exciton interaction and beyond the Markovian approximation 
for dephasing. Effects such as the polarization-dependent response of the excitons |I7|j20 



and signatures of bound biexcitons f2lH23f lead to more refined theoretical investigations 
beyond the mean-field approximation [pl)|, |I7| , [24] -|32f . 

Polarization-mixing in the nonlinear optical response is caused by excitation of two ex- 
citons of opposite spins which first occur in second order of the Coulomb interaction p? 



The mean-field description which treats the interaction between excitons to first order in the 
Coulomb interaction can account for neither the polarization mixing nor the formation of 



bound biexcitons. This point was clearly demonstrated by Combescot and Combescot [33 



who had stressed the importance of bound and unbound biexciton states for the excitonic 



ac-Stark shift as well as polarization effects [34]. The formation of para-biexcitons with 



singlet spin-states for the electrons and holes is one aspect of polarization-mixing which is 
most likely to be dominant for near resonant excitations of the fundamental exciton states. 
Moreover, even in the absence of bound biexcitons, the correlation in the continuum of 
two-exciton scattering states is also important [F2DJ,^D|. 

In studying the optical processes in semiconductor systems, while much can be learned 
from ensembles of non-interacting atomic transitions |35|-|37| interacting with the radiation 
field or interacting localized (dense) two- level systems [38] , the strong interaction and close 
proximity of the electrons in semiconductors provides a distinct avenue for new physics. It 
is well known that many-body effects lead to a renormalization of the external field (local- 
field effects) as well as to a renormalization of the interband transition energies (self-energy 
effects), depending sensitively on the density and dynamics of the surrounding electron-hole 
pairs P91H01. In the linear response regime, these effects assume quantitative importance 
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41] . However, the nonlinear properties of semiconductors in the weak nonlinear regime and 
in the high-excitation regime constitute challenging problem, where the electron interaction 
physics must be added to carrier non-equilibrium and physics of quantum optics. 

In the low-density or weakly nonlinear regime, i.e., to third-order in the external field, 
the dynamics of the semiconductor for near resonant excitation of the fundamental exciton 
resonances can be formulated in terms of a set of effective dynamical equations for the 
exciton polarization with nonlinear exciton-exciton interaction and space filling effects, which 
have been derived from the semiconductor Bloch equations (SBE) [p,4~2[,l3"l. These effective 



equations provide a useful and physical picture of the origin of the nonlinearities and the 
observed phenomena, for example in the theory of the four- wave-mixing ||, photon echo 
44| or the Rabi-oscillations in semiconductors [fl5|| . Being within a mean- field description, 



the effective equations can provide no more information than the full semiconductor Bloch 
equations in the mean-field approximation. 

The inclusion of biexcitonic effects as well as the possibility of polarization-mixing was 
first discussed in terms of phenomenological few-level models |16],[17|,|24| to explain the tem- 
poral dependence of the FWM-signal signal including oscillations as a beating phenomena 
between bound and unbound biexciton states. In more microscopic theories, it was first 
shown in the equation-of-motion method by Axt and Stahl E3] and later in a diagram- 
matic approach by Maialle and Sham |^7| that the semiconductor Bloch equations form a 
closed set of equations for the density matrix elements for any given order of the external 
field, depending on the initial state, which is usually the vacuum state of the electron-hole 
pairs. The so-called dynamics-controlled truncation scheme [^5| provides a starting point 
for a microscopic theory of the polarization effects and has been applied to up to fifth-order 



processes [23]. The inclusion of biexcitonic effects for actual applications is mostly treated 



in a restriction to bound-state contributions and a single two-exciton state contribution 



or perturbation theory in the Coulomb interaction f27|]28[ . The solution for the third-order 
susceptibility is inextricably bound to the solution of the four-particle problem. 

In this paper, we give a microscopic theory of exciton interaction effects in nonlinear 
optical processes based on the Coulomb interaction between electrons. The detailed account 



provides the derivation behind our results published earlier p0|,|4"6|| . The theory recovers the 



established mean-field results in the literature and formulates the rest of the interaction ef- 
fects, termed correlation, in a concise manner. Particularly striking is the resulting equation 
of motion for the third-order nonlinear optical response shown to be driven by a number of 
terms with clearly identified physical origins: (1) phase-space-filling corrections, which are 
due to the Pauli-blocking of electrons and holes, (2) exciton-exciton mean-field interaction, 
and (3) the correlation term which is expressed as a two-exciton force-force correlation func- 
tion. The derivation of the general equations of motion for any given order of the external 
field is given in Sec. J| in terms of the Hubbard operators, using a complete basis set of the N 
eh-pair states. The Axt-Stahl theorem |25| manifests itself as, for example, the Hilbert space 
of one and two electron-hole pairs being sufficient for the third-order nonlinear response of 
the semiconductor. Details of the commutation algebra are relegated to Appendix [A|. An 
alternative derivation in terms of the density matrix is not recorded here to keep the length 
of the paper within bounds. The part of dephasing which is due to the electron interaction 
effects is included in our correlation function and the rest of the dephasing due to other 
causes is treated phenomenologically. 
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The correlation function approach gives a unified description of all correlation effects. It 
naturally encompasses the recently observed polarization mixing and bound-state biexcitonic 
molecules. The exact two-exciton correlation treats these effects on the same footing as the 
two-exciton scattering states which will be shown to be equally important. It contains 
the exact spin-dependent Coulomb correlation among the four particles and determines the 
spectral weight of the biexciton states for the source-term of the nonlinear response. These 
properties can be demonstrated by a numerical example of a one- dimensional model system 
with the advantage of not making decoupling approximations of the correlation effects. (See 
Sec. PIT] ). Details of the model are described in Appendix [FJ. 

Section [IV] gives an example of the correlation effects on a nonlinear optical process: a 
three-pulse four- wave- mixing experiment. Extensive new results of the numerical evaluation 
of the simple model are qualitatively compared with experimental results. Exact numerical 
evaluation of the correlation function is confined to simple models. For more realistic models 
to the semiconductor systems, we need reasonable approximations. A number of these are 
investigated in section [Vj, including a microscopic expression for the excitation induced 
dephasing [[IT] and a brief comparison with the Boltzmann corrections in the quantum- 
kinetic equations |[48|| . We conclude with a summary of our theory and with a brief outline 
of future applications in section [VI|. 



II. AN EQUATION OF MOTION FOR THE THIRD-ORDER RESPONSE 



We take as the fundamental approximation that, in the absence of the light-matter 
interaction, the Hilbert space of the semiconductor model consists of disconnected subspaces, 
which can be labeled according to the number of electron-hole (eh) pairs in the many-particle 
states. Let |0) denote the trivial ground state with no eh-pairs present with energy co>o,o = 0. 
The one eh-pair subspace is the exciton subspace with states \E^ a ) with the quantum 
number n, a polarization index a and energy u>i )n)(T . Both bound and scattering states are 
included in n. The polarization index labels a specific transition which, if optically active, 
corresponds to the helicity of the light required to excite the eh-pair-states that form the 
exciton. For example, in Zincblende structures, four p-type valence band states with total 
angular momentum 3/2 are connected via an optical dipole transition to an s-type conduction 
band with spin degeneracy. Due to selection rules, the m = 3/2 (m = —3/2) electrons in the 
heavy-hole band are coupled via an optical transition with — (+) polarized photons to the 
s = l/2(s = —1/2) spin states in the conduction band. The m = l/2(m = —1/2) electrons 
in the light-hole band are coupled via an optical transition with — (+) polarized photons to 
the s = — l/2(s = 1/2) spin states in the conduction band. The spin-orbit interaction usually 
splits off a valence-band with total angular-momentum 1/2 and is neglected throughout this 
investigation. 

The next relevant subspace is the biexciton Hilbert space with a complete set \E$) 
of bound and unbound states. Here, we introduce a single index m to label the set of 
quantum numbers for the states with energy o^m- Even though not all the biexciton states 
are computed due to the many-body nature of the problem, we keep all the states as long 
as possible because occasions arise that the biexciton states as intermediate states can be 
re-summed by virtue of the completeness theorem, similar to the treatment of the ac-Stark 
shift [33]. Such a step would be lost in a common approximation which restricts from the 
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start to one or two biexciton states. 

The use of the subspaces of different exciton numbers as disconnected is implicit in pre- 



vious works p5| , p7| . The disconnectedness is an approximation because states with different 
exciton numbers can be connected by the Coulomb interaction. The most important conse- 
quence is the neglect of the electron-hole pair fluctuations which affects the ground state and 
the dielectric screening of the Coulomb interaction. In other words, we define the ground 
state of the semiconductor as a vacuum state with respect to the exciton annihilation 

B n>ir \0) = 0, (2.1) 

which means that no electron-hole pairs are present in the semiconductor ground state. The 
dielectric screening is approximately accounted for by the static dielectric constant of the 
semiconductor. 

We define a total a-polarization connected to an optical transition, which can be of 
arbitrary helicity, depending on the electronic states involved, 

P ff = /4£^ (2.2) 

k 

where the operator a creates a zero total momentum eh-pair with electron wave- vector k, 
hole wave-vector — k and polarization index a and \x a is the dipole matrix element between 
the electron and hole states, assumed to be independent of k. Completeness of the operators 
leads to an equivalent expression of the a-polarization in terms of exciton creation operators: 

P a = fJL* E «n,a5n l( T , (2.3) 
n 

where 

a n . a = v / l 7 $ n , CT (x = 0), (2.4) 

in terms of the exciton wave function at zero relative distance. The operator 

k 

creates an exciton state B\ a \fy = \E^' a ) with zero total momentum, energy iO\^ a and 
relative wave function 0k,n,CT in terms of momentum k. The combination of the electron 
band Ai and the hole band A2 determines the polarization index o = er(Ai,A2). The laser 
central frequency uj p is implicitly subtracted from the exciton energy when we transfer to 
the rotating frame. The biexciton energy uj2, m then contains a reduction — 2u p . Note that 
the factor a n ^ depends on the sample volume V of the system for bound exciton states 
and is nonzero only for exciton states with s-wave symmetry. We take care of this volume 
dependence and show clearly how the final result is indeed volume independent. 
Using the Dirac notation, we introduce the following Hubbard operators 

Xn,ol;M,i3 = |-EjV,a) \Em,p\ (2.6) 
which can be used, in combination with the completeness relation 
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1 = X] -^N,a;N,a , (2.7) 
AT, a 

to express the exciton operator as 

(2.8) 

N>l,a,/3 

The interaction of the semiconductor with a classical external laser field with central fre- 
quency Up and field-strength E(t) = J2a E a (t)e~ lu3pt e a + c.c. is given in the usual rotating 
wave approximation j35[ by 

^i = -E(^(*)^,a + h.c.) (2.9) 

with 

E n ^t) = fi a a* E a (t) (2.10) 



For comparison with previous work [43[ , this expression is the time-dependent renormalized 
Rabi frequency of a given polarization a and transition n (h = 1). The Hamiltonian of the 
semiconductor, from the disconnectedness of the subspaces, is 

(2.11) 

N,a 

which is equivalent to a multi-band microscopic Hamiltonian in second quantization. From 
the form of the interaction Hj it follows that the expectation values (X 0; 7v iQ )t can be ex- 
pressed as a power series in the external field 

(X0;N,a)t = E X ( N N : 2m \t) + 0{E N+ ^ +2 ) . (2.12) 
m=0 

The expectation value of a zero to iV-pair transition is at least of order N in the external 
field. This theorem has already been proven by Axt and Stahl [25|. An important relation 
can be derived by the identity for an arbitrary state | </>(£)) 

W)) 

= m)\E N , a ) (oW)) (<t>{t) |o> (E M ,pW)) (Xov); 1 . 

With Xq (t) = 1 from the initial condition of the semiconductor in its ground state |0), we 
find for the general expectation values 

{XN,a;M,P)t = (Xo;N,a)t (^0;M,/?)i (^0;C))t ■ (2.13) 

In order to calculate the a-polarization we consider the equation of motion for (B n ^ a ) t . 
Using the Hubbard operators it reads 

£(BnJt = Kn,. - iT)(B n , a ) t + EE ^ ;a ,p(X N , a ., N+1 , p ) t (2.14) 

01 N>1 a,/3 

- E n ^(t) + E ^2 E i N J;a,p(XN,a;N,p)t 
N>1 a,/3 
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with 

C nJ;a,l3 = (^N+1,/3 ~ ^N,a ~ ^l,n,a){E N ^ a \B ntCT \E N+ i J) (2.15) 

S nSa,/9 = H ^E^WwAEnA 1 ~ ™l,q,<x + nz^\E Nt p). (2.16) 
q 

We have introduced the dephasing due to degrees of freedom not included explicitly (e.g. 
phonons) in a phenomenological way with the effective parameter V. Using Eqs. (|2.12|) 



and (|2.13|) we see that Eq. Q2.14Q can be considered as a linear differential equation with 
a (trivial) first-order source and nontrivial source terms of third and higher order. In the 
following we restrict ourselves to the contributions up to third order. Then using Eqs. (p. 12 ) 
and fl2.13|) we see that only X^\ a and X2I have to be determined. As x[^l a (t) obeys 



Eq. ( |2.14 ) without the terms involving the summations one obtains 

= * f e-^^-^-^E^dt' . (2.17) 

J — oo 

The equation of motion for X^p reads 



i^X^ = (cu 2>p -iT xx )X^ (2.18) 
E E n " ta "( t )(E2,p\ B i" iCT "\Eiy,a'}x[]l l y 



with the biexciton phenomenological dephasing constant T xx . In order to write the second 
term on the rhs of Eq. (|2.14j) in a compact form we use the explicit result Eq. (|2.18|) for 
X\ p(t) in order to perform the summation over the biexciton quantum numbers (3. 

-i(ui2 t 0—iT X x)(t—t') 



n,a ,n ,a 

x E n ^ a n(t'){E 2 , p \Bl„ a „\E w )X[% >(Tl (t') (2.19) 



Using the identity 



\dt ( E /^ iCT »(*)^?^.'(*)e< Wl -''-'' +Wl -' CT '- 2ir )*J (2-20) 

\n',a';n ,cr / 



we obtain after a partial integration 

.(1) v(2), 



E<U,/^W (2-21) 

P 

o E ] {El,n,a\Bn,a {EI — UJ\^ h(J — ^l,n,5-) B\, j, \ E\^ n i a ) 

i i " " 
n ,a ,n ,a \ 
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p 

X 



xx^(t)x;;;, v ,(t) 

£j J— OO \ 

x x;J lCT/ (^0x;J /)CT 4^0e~ ^ ^ 1 •"'' CT '' + " 1 •"' CT '" 2^^)(^ " ^ ' ) ^ , 

We can also perform the /3-summation in the second term on the rhs of Eq. (|2.21| ) . For the 
derivative, we use the following identity, which holds for [B^z, B n>a \ = and H\0) = 

9 T ( (Ei t fi,&\Bn,<T {H — U)x, n ,a — ^l,n,a) e~ %Hl ' BK |-&L n ' iCT ) (2.22) 



X g*l a; l,n",<T"+ aJ l,n',<7'J T 



= -i(0|^ ;rij(T (T)^, iCT , n „ )CT „|0) e^WV^W,^ . 
Here we have introduced the "force" operator 

Dn,a-n,a = [Bn,a, {B n ,a, H}} (2.23) 

and the usual time dependence in the Heisenberg picture is given by D{r) = e lHT De~ lHl ~ . 
This allows to write Eq. ( p. 21 ) in a compact form defining a memory kernel 



<g£ (r) := (0\Dnfiw(T)Dl^„ y |0) , (2.24) 

for the second term on the rhs of Eq. ( |2.21|) . The matrix element in the mean-field contri- 
bution (first term on the rhs of Eq. (|2.21|) ) can also be simplified, which gives 

(El,n,&\B n ^ (H — iO\^ a — Wl,n,er) = \E\jh,a\ ([B n ,a, H] — COl, n ,a) (2.25) 

and we finally arrive at 



-if e~^F^-{t-f) 

J -co 



Here it is necessary to assume the relation T xx = 2T to use the identity Eq. ( |2.22[ ). The 
memory function in Eq. (|2.24|) is a four-point correlation function in terms of electron (hole) 
operators. As the operators A„ )(J (q) can be expressed in terms of finite center-of-mass 



S 



exciton operators B n>(T (q), F(t) can be considered a two-exciton correlation function. From 
the double commutator definition of the D operators, Eq. ( p.24j ) is a force-force correlation 



function fE) |. An explicit form of the operator D is derived in Appendix A using a specific 
semiconductor model. 

The first expression on the rhs of Eq. (|2.26| ) describes the correlations between excitons 
as in the usual mean- field semiconductor Bloch equations (MFSBE) It is only 

nonzero for excitons with zero center-of-mass momentum and with identical polarization, 
i.e., each pair of charged carriers (electrons and holes) must belong to the same bands, and 
consequently does not produce polarization mixing. For comparison with previous work on 



the MFSBE and earlier work on exciton-exciton interaction p0|-p2|, we introduce the matrix 
elements 

= (0|IV;^4v'4v' 1°)' ( 2 - 27 ) 

i i " " 



lUZ '* = (0\Dn,a;n,*Dl, ^ |0> . (2.28) 



As -Dn,<5-;n,o-(T)|0) = 0, the correlation function F(t) can be written as a time-ordered product 
and standard Feynman diagrams can be used e.g. to set up approximation schemes. From 
a diagrammatic analysis to all orders the rigorous polarization selection rule a + a = a' + a" 
can easily be read off |]27|| . The fact that the third-order polarizability can be expressed 
in terms of this correlation function depending on a single time difference is due to the 
simplicity of the semiconductor ground state Eq. ( |2.1| ) approximated by the vacuum state 
of the bound and unbound excitons. 

An additional contribution to the third-order nonlinear response is given by the phase- 
space filling factor, which is due to the Pauli blocking of electrons. This term is assumed to 
play a minor role in the low-density regime of optical excitations of semiconductors, but we 
include this contribution here to preserve the exactness of our expression to third order in 
the exciting field. Keeping the N = 1 contribution in Eq. (|2.16|) we find with Eq. (|2.13| ) 



E^S^W,/^ = E C<^^E n/ ,At) (2.29) 

a,/3 



(*8U(*) 



Vx (1) „ „(t) 

J l,n ,cr v ' 



The phase-space-filling parameter depends on the explicit exciton wave-functions, c.f. Ap- 
pendix [A]: 



q 

x (J5i jfii? |l - ni iqiCT + n 2 ^ u \E hn „ . (2.30) 

The exact third-order nonlinear polarization Plfl{t) is given by the solution of the following 
linear differential-equation with a complete set of source terms 

(j^ + ™W + r) P$(f) = -i J2 (4%a(t)Y (2.31) 



-_..//// 

n,a,n' ,<t' ,n ,cr 



^yri ,cr' ;n" ,cr" it, /,\y(l) /.\ , ± an' ,<r' ;n" ,<r" / f \y(l) (.\ 

- I f e-^F^-{t - t>)xU,At')xU, a „(t')dt> 
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This equation expresses succinctly the physical origins of the source terms which drive the 
third-order polarization: the first term in the curly bracket being the phase-space filling, the 
second the Hartree-Fock or mean-field terms (of first order in the Coulomb interaction be- 
tween excitons), and the last the biexciton correlation. The source terms require a solution 
of the linear-response problem of Eq. (|2.17 ). The above derivation of the third-order nonlin- 
ear response can be extended to higher order in the external field. However, the higher-order 
correlation functions involved no longer have the simple structure of the third-order response 
force-force correlation function. 



III. EXCITON-EXCITON CORRELATIONS 

The correlation function formulation of the last section provides us with a powerful 
framework to compute the correlation effects of the two excitons based on treating the 
interaction among the two electrons and two holes on a equal basis. In this section, we 
start with the nonlinear optical processes which lead to the excitation of the excitons and 
discuss some general properties inferred from a study of a one-dimensional system, where the 
correlation function can be calculated numerically for relatively large system sizes without 
approximating the effects of the many-body interaction. We calculate the spectral-function 
for the ls-exciton contribution to the third-order optical response for a one-dimensional 
semiconductor model. This is done in the frequency representation with the use of the 
Lanczos algorithm ||53|| . The details of the one-dimensional semiconductor model in real- 



space are defined in Appendix |B[ 

The selection rule which connects the helicity of the the exciting light to the spins of the 
electron and hole gives rise to different types of two-exciton excitations depending on the 
energy level structures. For a semiconductor with Zincblende structure, where spin-orbit 
leads to a four-fold degenerate valence-band maximum, Fig. [I] illustrates the three types: 

Type-I In heterostructures, the strain or confinement can lead to a considerable heavy- and 
light-hole splitting Ahi~ih- This allows for a selective excitation of heavy- or light-hole 
exciton states. 

Type-II In bulk systems, circularly polarized light excites a combination of heavy- and light- 
hole excitons simultaneously. 

Type-Ill Polarization mixing can be induced by a linearly polarized light. A special situation 
arises, when conduction band states are common for both excited exciton species, in 
contrast to the type-II excitation. A three-pulse four-wave-mixing geometry can distin- 
guish this analog of the A-transition in quantum optics from the decoupled transitions. 

Correlation among the electrons and holes depends sensitively on their angular momenta 
(spins) and, in turn, influences the polarization dependence in the nonlinear optical response. 
We sort out for each combination of spins the possible types of excitations listed above: 

• Opposite-spin excitons: 

If the two excitonic transitions belong to different conduction and valence bands, i.e., 
no single-particle states are the same, e.g., the polarization mixing of simultaneous 
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excitations of the m = —3/2 valence band to conduction-band s = — 1/2 transition 
with a positive circular polarization of the light-field and the m = 3/2 valence band 
to conduction-band s = 1/2 states with a negative circular polarization of the light- 
field respectively, this excitation condition is type-I. In the absence of a significant 
band-splitting, e.g., in bulk systems, a positive circular excitation will excite light- 
hole excitons from valence band m = —1/2 to conduction-band s — 1/2 states as 
well as from valence band m = —3/2 to conduction-band s = —1/2. This situation 
belongs to the type-II excitation. Although both excitons have equal-spin, the exciton 
energies and transition strengths are different and can be distinguished. In both cases 
above, (one type-I and one type-II), the hole and the electron states belong to different 
bands. The correlation among the four charged carriers is now entirely determined by 
Coulomb effects since no Pauli-blocking is encountered in the band. 

• Parallel-spin excitons: 

When the two exciton transitions have common hole and electron bands, their spins are 
parallel. For a type-I excitation, with only, say, positive circularly polarized excitation, 
only spin +l-excitons are populated. A three-pulse experiment can distinguish between 
the equal-spin and the opposite-spin correlations in the response signal. This point 
is discussed in the context of the four-wave-mixing experiments in the next section. 
The most prominent feature of the equal-spin correlation is the absence of a bound 
biexcitonic molecule. This shows up clearly in the one-dimensional model. A type-II 
excitation will mix the opposite- and parallel-spin correlations and the analysis of the 
FWM geometry is more complicated, due to the natural interference of the individual 
contributions. 

• Coupled-spin excitons: 

The type-Ill excitation arises when either the electron or hole states but not both 
share a common band. For example, linear excitation of heavy- and light-hole exci- 
tons in GaAs will induce transitions from the m = —3/2 valence states with photons 
of positive helicity as well as transitions from the m — 1/2 valence states with photons 
of negative helicity to common conduction-band states with s = —1/2. The corre- 
sponding correlation function is similar to the parallel-spin case in the sense that the 
conduction-band electrons are subject to the Pauli exclusion. With a finite heavy/light 
hole splitting, this type-Ill excitation is sensitive to the time delay between the two 
exciting pulses in a three-pulse FWM experiment and can, thus, be distinguished from 
the parallel-spin type-I case. 

Fig. ^] shows the correlation function spectrum for the diagonal ls-exciton contribution 
calculated in a one- dimensional semiconductor model. The opposite-spin case (solid line) 
and the parallel-spin case (dashed line) are shown. The zero of energy corresponds to the 
energy of two non-interacting ls-excitons, which is spin independent. Resonances at negative 
energies to correspond to bound excitonic molecules (biexcitons). In this case, the binding 
energy of the biexciton is approximately 1.5 meV. In Fig. |2| (a), the masses of electrons and 
hole are identical (positronium-limit). The bound-state of the para-biexciton is the most 
significant feature at low energies, since the binding energy is expected to be much smaller 
than the usual excitonic binding energies, which is 10 meV in this model for the ls-state. 
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The parallel-spin case has more pronounced spectral weight at lower positive energies, but a 
bound-state is not expected to exist. The spectra have a maxima at higher energies before 
dropping to zero. The bandwidth of a single-electron band is 50 meV. It is surprising that 
the spectrum is almost zero above the free electron-hole pair-state bandwidth of 100 meV. 
However, a high energy-resonance is visible. 

This distant resonance moves to lower energies in Fig. § (b), when the ratio of the 
electron-to-hole mass is reduced to m e /mh = 0.15, found in semiconductors like in GaAs. 
The reduced mass of electrons and holes is kept constant. More spectral features appear 
in the lower energy regime. In principle, more bound-states should appear for negative 
energies. A dip-like structure can be found at almost the same position for the opposite- 
spin and parallel-spin case. In the case of opposite-spin correlations, we find increased 
spectral weight for smaller positive energies. This can be easily understood, if one eh-pair 
is in the ls-exciton state, whereas the second eh-pair is quasi-free (dissociated). 

This can be seen in Fig. [2] (c) for infinite hole mass (molecular limit). The opposite- 
spin case has pronounced spectral features at low energies. More resonances appear, if the 
spectral broadening in decreased. The resonances can be classified as (1) bound excitonic 
molecules at negative energies, which simply are the ground-states of two electrons in a static 
potential of two heavy holes with varying distance, and (2) scattering resonances, which are 
anti-binding states with molecular character. These anti-binding states do contribute to the 
nonlinear optical response. For example, the first large resonance at positive energies results 
from the configuration where the two holes are extremely close together. This is similar 
to the configuration of a helium atom from the view of the fast moving electrons. The 
spectral weight of the force-force correlation function favors states with small distance of 
the charged carriers, as can be seen from the real-space representation of the state -D^O) in 
Eq. ( |B13| ). A third resonance is clearly visible at the same energy position for both spin cases 
at uj pd 10 meV. This resonance is the largest feature in the case of parallel-spin excitons. 
The analysis of small systems of up to 10-sites strongly suggests that this feature originates 
from Coulomb correlation of electrons and holes in the anti-bonding state, where the holes 
are located on neighboring sites. This explains the feature in both spectra. The position 
of the resonance is only weakly dependent on the on-site Coulomb interaction. We have 
verified that the spectral resonance from a state with a single ls-exciton and a dissociated 
eh-pair gives much a smaller contribution at this energy. The dip for parallel-spin excitons 
in Fig. ^ (c) goes almost to zero and separates a small band of biexciton scattering states 
with weakly repulsive interacting pairs of ls-excitons from the anti-bonding resonance. For 
opposite-spin excitons, a spectral hole emerges at zero energy. Spectral weight from the 
product state of two ls-excitons with zero center-of-mass momentum is recovered in the 
mean-field contribution to the nonlinear optical response. 

The force-force correlation spectrum has a large number of spectral features at the exact 
biexciton energies, which include bound-state as well as scattering-state contributions. The 
correlation function approach treats these states on an equal footing. We expect a similar 
behavior also for more realistic 2d- or 3<i-models. 
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IV. APPLICATION TO FOUR- WAVE-MIXING 



In this section, we investigate the correlation effects in four-wave-mixing (FWM) exper- 
iments with semiconductor heterostructures in the low-density excitation regime where the 
third-order theory in the exciting field is valid. Polarization mixing of electrons and holes 
with different spin can be induced and probed using cross-polarized laser excitation. With 
the proliferation of new experiments which aim at probing correlations, we use, as illus- 
trations of application of our theory, two specific experimental situations with the simple 
type-I excitation. The first paradigm experiment I^Ofl , by resonantly exciting excitons in 
quantum-well structures, clearly demonstrates the signature of polarization-mixing in these 
systems for the ultrafast nonlinear response, even in the absence of biexcitonic molecules. 
The actual numerical simulations are performed with the quasi one-dimensional semiconduc- 
tor correlation function (Sec. 3) to model the spin-dependent effects and are not intended 
to quantitatively describe the experiment. As a second example, we study the effects of 
bound and unbound two-exciton states and discuss the "beating" phenomena of the biex- 
citon resonance for type-I excitations, which was identified as quantum-beats |E| between 
bound and unbound biexciton states and which our calculation shows to be a ringing of the 
bound-state resonance alone. 

The typical experimental setup is sketched in Fig. [I] for a three-pulse four-wave-mixing 
geometry. This experiment leads to polarization mixing when pulses (2) and (3) have op- 
posite helicity. In this case, equal populations of excitons with opposite angular momenta 
are excited for type-I excitation. In the mean-field description, no third-order polarization 
p( 3 )(k/) exists to diffract probe pulse (1) in the kf = k3 + k2 — ki direction for any time 
delay T. However, beyond the mean-field theory, correlation effects between opposite-spin 
excitons cause polarization mixing. 

To allow for an analytical discussion, we consider the limit of ultrashort light pulses with 
identical central frequency, which simplifies most of the following calculations. In this case, 
the exciting laser field is given by 

E(t) = ^E(k,) ( 5(t + r J ). (4.1) 

j 

We emphasize that a delta-pulse approximation is used only in the time integration and 
not in the frequency integration since it gives a infinite broad spectral width. The laser 
pulse in the propagation direction k 3 interacts with the sample at time t = —Tj. The 
corresponding area of the light pulse is given by E(kj). The summation index j labels the 
different pulses involved the multi-wave experiment. The first-order polarization of a single 
exciton transition is easily calculated from Eq. ( p. 17 ) 



= E e-^-^ t+ ^Q(t + tAE^). (4.2) 

3 

Here, E a (kj) = E(kj) ■ e CT is the projection of the laser field in the propagation direction 
kj onto the polarization unit-vector e CT of the (n, a) transition. The dephasing T a may 
depend on the polarization of the transition. In the exciton picture, the total nonlinear 
polarization to third-order in the external fields is given by Eq. (|2.31|) . The contribution 



from the phase-space-filling, well documented in the literature |p4|,^|,|3T| , plays only a minor 
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role in the low-density limit. We focus on the remaining contribution from the mean- field 
part and the genuine correlation which can be treated on equal footing. Correlation lead to 
the following nonlinear complex polarization for the transition (n, a) 

P&V) = -^> n ,,e-^-^ £ (aS<P) (4.3) 



n,a ,11' ,a f ,n ff ,<r" 

ii,J2>i3 



{ e(t + r il )e 21 e 3 2C' 12 (t) + e(t + r i2 )e 12 e 32 c 22 (t) 
+ e(t + r h )e 31 e 23 c 13 (t) + e(t + r js )e 13 e 32 c 33 (t)} . 

with Qm = ©(Tfc — t~i) and an exciton wave-function dependent factor 

a = an,a< v < v », (4.4) 

with a n>a defined in Eq. (|2.4|) . The exciton- label dependence of a on the left is understood. 
The external field and helicity dependence is contained in the factor 

S = iiaE~(k j J/j la ,E a ,(\i h )iJ l(T r,E tT r,(k h ). (4.5) 

A general phase <fi is due to the delay between the short pulses 

(f) = e^^^+^^nQ-^ri' .a'-^a'^Q-^n'' ^''-^a' 1 )^ . (4.6) 

The nontrivial part of the polarization dynamics in contained in the time-dependent function 

(t > -T a ,r a < r b ) 

e i(«n, <r +wft,ff-a; n /, ff /-a; n //, (r //)f e -(r (r ,+r (r „)f / p(r)dTdt'. (4.7) 

-Ta JO 

Eq. ( f4.7|) contains the complete Coulomb correlation beyond the mean-field approximation. 
The integral kernel is the memory function, which can be calculated with the knowledge of 
the exciton-exciton correlation function, 



(r) = e-^'^"^F(r). (4.8) 

For the numerical calculation, we always use the resonant condition for the temporal evo- 
lution of F(t) in Eq. ( |4.8j ), i.e. uj n i + u n n = 0. For non-resonant excitation, the detuning- 
dependent integrand in Eq. (|4.7| ) has no phase dependence for the diagonal contributions of 
the correlation function. 

We now consider the response from the non-degenerate exciton ground-states (Is) only, 
with near-resonant excitation of the central laser frequency, i.e. ou n =is,a = and neglect 
the transition label n = Is. We also set T a = V for convenience. This simplifies the general 
expression Eq. ( J4.3| ) and Eq. ( |4.7| ) considerably: 

C a>b (t) = f e~ 2Tt ' f +T " F(r)dr (t > -r a , r a < r b ). (4.9) 

J-Ta JO 

The simplest experiment which is able to distinguish equal-spin from opposite-spin cor- 
relation is the three-beam experiment of Fig. |3], where the two pulses (2) and (3) interact 
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with the sample at the same time. Assuming that lasers (2) and (3) interact with the sample 
at r 2 = r 3 = we define the delay time t% = T. The correlation function is now diagonal 
with respect to the exciton indices and has the spectral representation with f m = (0\D\E m ) 

Hr)=Y,\U 2 e- luJmT (4-10) 

m 

where the m-summation includes all contributions from bound (u m < 0) and unbound 
{uj m > 0) biexciton states \E m ) and the implicit transition indices are understood. We find 
for Eq. (f4.9|) with r a = T, r b = for negative time delay T < and t > —T, 



„ J-T Jo 



t' 



e -(2r+iw m )i _ e (2F+iw m )T e ~2Tt _ g 2rT ' 



El/^l 2 • ( Ion or • ( 4 - n ) 



The second term of Eq. Q4.11p is simplified by a sum rule 

EI/JW = ^ ( 4 - 12 ) 

m 

which can easily be derived from the usual Lehmann representation of the correlation func- 
tion and using Eq.( 2.27| ). The parameter (3 is the mean-field exciton-exciton interaction 



parameter with four identical exciton indices (Is). For the parallel-polarized case, this 
contribution is canceled exactly by the explicit mean-field contribution to the third-order 
nonlinear polarization. This can be seen by inspection of the second and third terms on the 
rhs of Eq. (|2.31|) , when we set F(r) = i(35(r). For opposite-spin excitation, this parameter (3 



is zero. In the following we define C(t, T) to be the result of Eq. (4.11) after the cancellation: 



C(t,T) := E \f m \^^ T \ 1 (4.13) 
and, for comparison, the mean-field part 

t -2T{t+T) _ i ^ 

C MF (t,T) := -t(3e^ J £ i j . (4.14) 

The sum-rule Eq. ( [4.12|) does not imply that the mean-field instantaneous contribution in 
the nonlinear response has completely disappeared in Eq. ((4.13[). Mean- field and correlation 



contributions are treated here on an equal footing as part of the Coulomb interaction of 
the charged carriers. In the response function Eq. ( 4.12 ), the mean-field contribution is 



recovered in the large dephasing limit, i.e., T >> uj x . In the following, we assume equal 
field strength for all pulses with real amplitude E a . We evaluate the signal according to the 
following general rules: 

1. Spatial dependence. 



The external sources, Eq. (4J3), for the nonlinear polarization Eq. (4.3) have to be 
selected with the correct spatial phase dependence. For the signal in kj-direction, 
the index combinations k J2 = k 2 , k J3 = k 3 and kj 2 = k 3 , kj 3 = k 2 are possible with 
kj-j = k x fixed. 
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2. Time dependence. 

The set {ji, j 2 , J3} determines the type of response term and the corresponding 
0-functions. This will determine the temporal details of the signal, depending on the 
time-order of the incoming pulses. 

3. Helicity dependence. 

The polarization of the transitions and the helicity of the exciting fields determine the 
correct type of correlation function and field amplitude projection E a . 

4. Transition dependence. 

Perform the summation over the exciton quantum numbers after the polarization de- 
pendence is determined in the above steps. Symmetry arguments can be used to reduce 
the actual number of terms. 

Applying rules (1) and (2) in the case of a near-resonant excitation of the heavy- 
hole/light-hole ls-excitons, we find for the time-resolved nonlinear polarization, using 

^(*) = -k E e-^- ir )V^+ ir >V^(ki) (4.15) 

1 ~ it 

x ^E (T /(k 2 )^/£; (T //(k 3 ) +/v£ <7 /(k 3 )/v/£>(k 2 )) K (t, T ) 



with 



K (t, T) = Q(t)Q(T)C(t, 0) + Q(t + T)Q(-T)C{t, T). (4.16) 



The helicity of the diffracted polarization depends on the individual contributions, as mixed 
in Eq. fl4TTJ). 

We discuss in more detail a type-I excitation. Corrections due to phase-space filling are 
neglected. The three-pulse experiment in Fig. |3] can clearly distinguish between correlations 
between opposite-spin excitons and parallel-spin excitons. In the cross-polarized configura- 
tion, pulses (2) and (3) have opposite helicity. The response in kf-direction has opposite 
helicity with respect to pulse (1). This is a consequence of angular-momentum conservation. 
In the case of equal helicity, all three pulses and the response have identical helicity. In both 
cases, the main difference comes from the type of heavy-hole correlation function, which en- 
ters the calculation of C(t,T). Since the response comes from identical optical transitions, 
besides the helicity, we set \i a E a e {0, 1}. The responses in both polarization configurations 
differ mainly in the type of correlation function to be calculated. For the ±-heavy-hole exci- 
ton transition, we find for kx — > a, k 2 — > a and k 3 — > a, in the cross-polarized configuration 
and kj — ► ex in the co-polarized configuration (a = —a) 

pQ)u\ = _ ie -i(u>a-ir)(t+T) K r t ^ T \ J e 4( ^ +ajCT)T cross-circular, U-lT) 

1 e i i< jJ °+u<*) T co-circular. 

For resonant excitation, i.e., uj a = 0, the time- resolved phase of the polarization can be read 
off P( 3 ) = {P^le 1 * H, 
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$w ._ e(i + r)e(-r )ta ^(i^) (4,8) 

in the case where probe pulse (1) comes after the excitation with pulses (2) and (3), T < 0, 
and the case where pulse (1) precedes the excitation, i.e., T > 0. For small times t after 
the excitation, the instantaneous phase-space-filling term leads to a 7r/2-phase shift of the 
polarization with respect to the external field E, as discussed in |18 |. 
The observed quantities are the time-resolved (TR) intensity 

I^(t,T) = \P! r 3 \t)\ 2 

= e -2T(r+*)|«(t s T)| 2 (4.19) 

and the time-integrated (TI) intensity 

roc 

/(3)(T) = / \P®(t,T)\ 2 dt. (4.20) 



We have performed numerical simulations using a one- dimensional extended Hubbard model 
with long-range Coulomb interaction as defined in Appendix [FJ. We present in this section 
the exact numerical calculations for this simple model and consider in the next section various 
approximations involving truncation of the summation in Eq. (|4.13|) for more complicated 
models. For resonant excitation, the ultrafast polarization dynamics is strongly affected by 
the relation between the exciton dephasing parameter T and the Rydberg energy u x , since 
the detuning is zero. The Rydberg energy does not appear explicitly, but for the Coulomb 
interaction U ~ u x holds. The results can be compared with few-level models on FWM P| 



The source term C(t, T) of Eq. (|4.13|) plays a central role in the nonlinear response, since 
it determines the nontrivial polarization dynamics. Many-particle correlation leads to a 
dynamical structure which is absent in a simplified non-interacting two-level system. Fig. |] 
shows the typical source term for the parallel-spin case. The mean-field source term shows 
a finite rise time which corresponds to the finite rise time of the time-resolved nonlinear 
polarization signal, roughly the dephasing time T 2 . For larger times, the nearly constant 
source term leads to an exponential decay of the resulting TR-signal in Eq. ( 4.15 ). The mean- 



field picture is considerably changed when the exact correlations are taken into account. 
Fig. ^ shows the following characteristic features: (1) an increase in the rise time of the 
signal compared to the mean-field approximation, (2) the signal exhibits a phase dynamics, 
and (3) the asymptotic value is complex and differs considerably from the mean-field value. 
Only in the extremely large dephasing limit, not shown in the figure, the correlation result 
approaches the mean-field value if3/2Y. 

Fig. |5| shows the typical source term for the case of opposite-spin correlation. The 
existence of a bound-state biexciton has a strong influence on the nonlinear response. Oscil- 
lations with the biexciton binding frequency of the single bound-state in the one-dimensional 
model are visible in Fig. [5[ The energy denominator in Eq. ( 4.13 ) favors low-energy reso- 



nances, i.e., isolated bound-states {u m < 0) and the low-energy scattering-states continuum. 
It is important to note that the oscillations decay with twice the polarization decay-time 
T 2 = T _1 in the approximation of Sec. 0. In the TR-signal, the oscillating contribution to 
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the signal should, therefore, be fairly small on the decaying part of the signal. The second, 
more important, observation is that, from Eq. ( 4.13Q and the spectrum in Fig. f|, the biexci- 
ton resonance alone is responsible for the oscillations. This is a ringing phenomenon, which 
is different from the usual quantum-beat picture which is suggested by the few-level model 

mm- 

Fig. |6] shows the results for the time-resolved polarization for zero delay T = and 
weak dephasing. The solid line shows the mean-field calculation, which obviously gives a 
poor approximation for resonant excitation and small dephasing. The exact result for cross- 
polarized circular excitation shows the ringing of the intensity signal with the biexciton 
binding frequency, which is already present in the source term in Fig. [5| The signal for 
co-polarized excitation is of the same order of magnitude but shows no oscillatory behavior. 
The signal peaks at roughly T 2 /2, which is less than estimated previously p0| . 

In Fig. |7|, a shorter dephasing time of T 2 = 0.5 ps is assumed. The mean-field result 
looks better in comparison with the co-polarized signal. The exact signals have a decreased 
rise time of the maximum, which can be explained by an additional dephasing mechanism 
due to the superposition of the continuum of two-exciton states, which leads to a natural 
intrinsic decay, similar to the effect of an inhomogeneously broadened system. 

The time-integrated signal in Fig. ^ shows the effect of finite delays between the exciting 
pulses (2) and (3) and the probe pulse (1). The result of the numerical simulation is in very 



good agreement with the experimental results by Wang et al. ||20|| , who have performed a 
three-pulse FWM experiment on a GaAs quantum-well system. The figure shows the strong 
decrease of the co-polarized signal for negative time delay, which is due to the enhanced 
intrinsic dephasing of the continuum of two-exciton scattering states. The cross-polarized 
signal is stronger for negative time delay, which indicates the effect of quasi-bound excitonic 
molecules. The spectral weight of the correlation function in Fig. ^| (b) for parallel-spin 
excitons is more enhanced in the low-frequency regime. This leads to a stronger total signal 
in Eq. ( |4.15| ), which can be reproduced for different dephasing parameters. 

Fig. |^ shows the time-integrated intensity of the FWM-signal for co-polarized circular 
excitation for different dephasing times. For short dephasing time, we recover the well-known 
mean-field behavior for homogeneously broadened systems, which predicts a rise of the signal 
~ T 2 /4 and a decay of the signal ~ T 2 /2 This can be explained by simply counting 

the number of polarization-waves which are present before the nonlinear signal is emitted. 
For positive time-delay, even for a longer dephasing time, the decay ~ T 2 can be observed 
because correlation effects influence only the strength of the time-integrated signal. For 
negative time-delay, significant deviation from the ~ T 2 /4-law is found. The probe pulse 
interacts with the delay —T > and fast decaying modes of the correlations cannot be 
sustained. The calculated results show a smooth transition from a steep rise near zero delay 
time, where correlations with the fast modes of the spectrum of two-exciton scattering states 
are important, to a regime where low- energy modes dominate the response. The latter again 
shows the asymptotic mean- field like ~ T 2 /4 dependence. 

The most prominent feature of the cross-polarized response in Fig. ID is the modulation 
of the signal with the binding-frequency of the bound-state molecule at negative time delay. 
This biexcitonic-effect has been observed experimentally by various groups |16],|22| and has 
sparked much theoretical effort to improve the mean-field theory of the semiconductor Bloch 
equations. This signature clearly shows the importance of correlations, which cannot be 
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neglected for the resonant excitation of the ls-exciton in semiconductor heterostructures. 
No signal in the cross-polarized configuration is predicted by the mean-field approximation. 
For shorter dephasing times, the modulations disappear very quickly and the signal shows 
similar behavior compared to the co-polarized geometry. We also observe a decrease of the 
modulations, if the pulse-width is increased. The biexcitonic-modulations (if a bound-state 
molecule exists) can only be observed for sufficient short laser pulses and weak dephasing 
of the coherence (polarization) in semiconductors. However, none of the above conditions 
is necessary to observe correlations due to the two-exciton scattering-continuum, which 
is always present and gives the main contribution, if bound-states are absent because of 
impurity-scattering, interface effects, etc. The correlation function approach gives a unified 
description of the observable effects. In principle the same treatment can be applied to type- 
II and type-Ill excitations, but the large number of terms involved renders the computation 
quite tedious. Much work still has to be done. A discussion of further specific nonlinear 
optical experiments where correlations are involved is in progress. 



V. APPROXIMATIONS AND COMPARISON WITH RELATED APPROACHES 

In considering possible applications of the correlation function approach for an improved 
treatment of the dynamical nonlinear response, we need a tractable response theory which 
takes into account, for example, (i) non-resonant, above band-gap excitations where electron- 
electron scattering becomes important and where the Markovian approximation is no longer 
valid on short time scales, (ii) coherent and non-coherent scatterings with LO-phonons, and 
(iii) applied magnetic field in heterostructures. Some theoretical work has already been done 
in deriving scattering-rate corrections with memory kernels for the SBE [[56] , in LO-phonon 
corrections pTLpSl], and in high magnetic fields ||32||. In semiconductor heterostructures 



or bulk systems, the calculation of the force-force correlation function, even for resonant 
excitation, is an enormous numerical task because of the four-body problem involved. For 
systems beyond the simple models for which the exact calculations are possible as described 
in the last two sections, we develop various approximation schemes to incorporate correlation 
beyond the mean-field level in the dynamical optical response for low-density excitation. 



A. Excitation induced dephasing (EID) 

EID corrections for the SBE have been discussed by Wang et al. [fTJ] in an application to 
FWM, where a phenomenological, density dependent and k-diagonal dephasing-parameter 
was introduced. We can derive a similar correction from the exact third-order contributions 
to the nonlinear response. In this case, the contribution originates entirely from exciton- 
exciton correlation. 

Extracting the fast time-dependence of the linear polarization in Eq. ( |2.17| ) 

*8U*) = e-*" 1 -" '*SU*) (5-1) 
and the correlation function, c.f. Eq. ( f4.8|) , we find for the correlation part of Eq. (|2.31 ) 
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f e-^>F%&"* (t - f (f a (t')dt' (5.2) 

J — oo ' ' 

J —oo ' ' 

~^ ln,a;n,a l,n' ,a' \ L ) ^ i >n " j(7 " \° J ■ 

Hence, we have identified: 

. . II II rOO „ i i " *> 

~n,a;n ,a _ / -2Vr pn ,a ;n ,a ( \J (r o\ 

ln,a;n,a ~ e r h,a;n,a \T )(1T , 



a microscopic expression for the EID for the low-density optical regime. The imaginary 
part of 7 renormalizes the mean-field parameter f3. The real part of 7 leads to a dephas- 
ing of the nonlinear polarization, which can be seen from Eq. ( |2.31| ). In addition to the 
phenomenological treatment ||47|| , polarization- mixing is now also automatically take into 



account. 



B. Short-time memory approximation 

In the limit of large dephasing, a short-time approximation of the memory kernel in 
Eq. is valid by Eq. ( ggg ), 



tit! . . !l II 

in' ,0" ;n ,cr / \ n 'o :n ,a 



" (r) -> 7&X ' + 0(r) • (5.4) 

Thus, we obtain a correlation-modified complex mean-field parameter, which exhibits po- 
larization mixing and leads to EID: 



n ,ct ;n ,<r 

~n',a';n ,a" lh,a\n,a /r r\ 

We observe that, for diagonal contributions, e.g., n = n' = a and n = n" = b, the expression 
Eq. ( |5.5| ) is positive. Hence, the nonlinear polarization for, say a, has an effective dephasing 
which depends on the density of species b given by 

r a - r a + -^-x;x b > o. (5.6) 

J- a + 1 6 

Explicit results for the expectation value 7^ are given in Appendix [A]. 



C. Non- interacting excitons approximation 

We propose an approximation scheme for the correlation function F, where we replace 
the time-evolution for the D-operator Eq. ( |B13|) with the full Hamiltonian H of the biexciton 
subspace (N p = 2) with the free time-evolution of excitons. For more explicit results, we 
use the D-operator representation of Eq. (|A13| ): 
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Fn£n^ ,a ( r ) ^ q ^-q'a„, a (q)a fiia /(-q)a n / i/ 3(-q / )*a n //^(q / )" 



x (OlS^CeOS^C-qjSj^C-qOB^CcjOlO) (5.7) 

x e -*( w / 3, t T'(q')+^ /3 ', ( T"(q'))^ i 



This expression can be simplified further. We note that the expression Eq. (577) is a non- 
perturbative result which is exact to second order in the Coulomb interaction but properly 
includes all Coulomb effects to form excitons. If no bound-state biexcitons are present, this 
expression should give a reasonable approximation beyond the short-time memory approx- 
imation. The commutator in Eq. (|5.7Q is calculated in Eq. ( |A14j) . Since we are interested 
in the slowly varying contributions to F(r), we restrict the remaining exciton summations 
in Eq. ( |5.7| ) to diagonal contributions, i.e., (3 = n" and f3' = n'. The summation cannot 
be performed analytically for finite r. According to Eq. (|A14|) , we can split the correlation 
function into three contributions. The bosonic part gives 

p(l) _|_ _p(2) _ ^2 e -i(a Vv (q)+c iV , (7 ,,(q))T ^ 

q^O 

+ ffi,ri (Q) fn,n" ( _ l) fn",n" ( — *l) * fn' ,ri ( 0) * ^cr' ^cr,c 



with the functions / defined in Eq.(A12). Due to the composite character of excitons, the 
third contribution is 

F (3) = _ J- f7 q f/ k _ k , e - i( ^'.-' (q)+ ^".-" (q))r (5.9) 

q,k,k' 

X \°o-,cr" 5-,cr' o-,cr' S,cr"in",n"l VU Jn',n'VVU 

+ 4«(-q)/Jn'(q)C4i'ftC") 

X ( $ f t ,k'+q,a - 3>ft,k+q,CT)($n,k,<T _ *R,kv) $ n",k'+^„ K(T „q, CT " < I ) n ',k+^,q, CT ' • 

The Kronecker-delta is restricted to the conduction (c) or valence-band (v), respectively. 
Eq. ( |5.7|) is a non-perturbative result and corresponds to the summation of a class of di- 
agrams in the expansion of the correlation function. This can be useful in calculating the 
response from higher dimensional semiconductor models, which will be addressed in future 
work. The low-frequency range of the spectrum is approximated by Eq. Q5.7| ), which corre- 
sponds to the long-time behavior of the correlation function and is dominant in the nonlinear 
optical response if no bound-state molecules are present. Higher-frequency modes naturally 
decay very fast, which affects the ultra short-time memory. Eqs. ( |5.S| , |5.9|) can be calculated 
numerically with the exciton wave-functions as input. 



D. Comparison with collision terms in the Boltzmann approximation 

Finally, we want to approximate the correlation to second-order in the Coulomb inter- 
action. This will enable us to make a qualitative comparison of the correlation-part of the 
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nonlinear optical response with semi-classical Boltzmann-Equation approaches. This is done 
here to third-order in the external field. We expand the polarization for a transition a and 
Bloch vector k in terms of excitons and inspect the source terms which gives the correlation 
corrections only. Since the force-force correlation function is already second-order in £/ q , we 
can use the non-interacting linear polarization Eq. ( j2.17| ) 

X^ a (t) = w£ f (E^\e^ H ^W\V,<r)E {l!)dt! (5.10) 
k J -°° 

k J -°° 

k 

The correlation source term in Eq. ( |2.31| ) to second order in the Coulomb interaction is 

:=E^k,.Pn, CT (t) (5.11) 

E ®n^nM,°K>te,°'®*n"te,°" (PSA*)] 



scatt 

n 
1 

~ 2 



- ... a a 

n,n,a ,n ,(7 ,n ,a 
k l> k 2. k 3 



x f e- 2r « 1 - 1 ')^ / ( < -i') P « ff ,(t')^,,(0^'. 

•J — oo J ' 

The exciton summations can be performed to yield 



k l> k 2> k 3 



= E 

ki,k2,k3,q,q' 

x (0\cl uai c kl+ ^A,^ + n^ 

x (^(*)*p£^(Op£V -pff-^WpS^co^ 
+ p k ?,^)>ki q ><(o^ q> » + pJUritw^/ 

-(k-k-q), 

with the non-interacting Hamiltonian if being used for the time evolution. The matrix 



element in Eq. ( 5.12 ) is divided into 4 contributions with a total of 64 terms, which can 
be simplified by symmetry arguments. Note that the bosonic contributions to the matrix 
element yield a 5 q ,±q' factor. This reproduces the third-order limit of the collision terms 
which were previously derived ||56|| . This also corresponds to the two-loop diagrams in 
the diagrammatic approach ||27|| . However, the true fermionic contributions also appear, 
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which are relevant for the parallel-spin and coupled-spin case and correspond to the one- 
loop diagrams. These terms might not be present in the usual semi-classical treatment of 
the scattering rates. The correlation function approach naturally incorporates these effects 
on the four-particle level and gives the exact low-density results. 



VI. CONCLUSION 

In this paper, we have presented a unified theory of exciton-exciton interaction effects in 
the third-order nonlinear optical response, using a correlation function approach J3(|. The 



electronic problem (dynamics of the four interacting particles) is separated from the non- 
linear optical problem. Furthermore, the correlation effects beyond the mean-field terms is 
explicitly represented by a two-exciton force-force correlation function. By means of this for- 
malism, we are able to investigate the role of exciton-exciton correlations in the third-order 
polarization in an application to resonantly excited heavy-hole excitons in a semiconductor 
quantum well. The correlation functions are calculated numerically for a one-dimensional 
semiconductor model with long-range Coulomb interaction, without perturbative approxi- 
mation. Their spectra exhibit isolated resonances due to bound-state biexcitons and con- 
tinuum of two-exciton scattering states. Additional, more pronounced features appear for 
decreasing mass-ratio of electron to hole. 

A three-pulse FWM configuration can distinguish between parallel-spin and opposite-spin 
correlations. For co-polarized excitations, we find a significant deviation from the expected 
mean-field ~ T2/4 rise of the time- integrated FWM-signal for negative time delays between 
the pulses in the weak dephasing regime. This can be explained by an ultrafast intrinsic 
decay of correlations due to two-exciton scattering states, since bound-state molecules are 
absent in the case of parallel-spin correlations. For strong dephasing, we recover the well 
known mean-field results 0. The signal in the cross-polarized configuration is dominated by 
a modulation with the binding frequency of the bound-state biexciton in the system, which 



has been observed in various experiments [16.22]. These oscillations appear as a ringing of 



the biexciton mode in the time-resolved signal, as distinct from true quantum-beats. For 
positive time delays, the third-order response is close to the mean-field T 2 /2-decay behavior. 
Polarization- mixing for cross-polarized excitations is purely a correlation effect ||27|| , absent 
in the conventional mean-field approach using the semiconductor Bloch equations PJTl|]. 

We have also derived generalized effective equations of motion for the polarization for 
laser excitation near the fundamental exciton resonance or the exciton continuum states. 
This is intended for a future application of the theory to more realistic semiconductor 
quantum-wells. The exact third-order equations are compared with the results of Boltzmann- 
type scattering corrections to the quantum-kinetic equations. The effective parameters en- 
tering the dynamics can be calculated in second order of the exciton-exciton interaction, 
provided only that the exciton states are known, which is possible for a large number of 
semiconductor systems. 

The correlation function approach can be generalized to include additional interactions, 



such as spin-flipping scattering processes or the coupling to LO-phonons |p8| . A further 
application of the correlation theory to spin-beating phenomena in diluted magnetic semi- 
conductors for the pump-and-probe configuration |59| has been made [60 |. 
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APPENDIX A: COMMUTATOR ALGEBRA FOR THE FORCE OPERATOR 

In this Appendix we derive the operators and parameters defined in section |T| in terms 
of a semiconductor Hamiltonian, H = Hq + U, with an independent-electron part: 

H = J2 e k, s c kjS c k , s , (Al) 

where c k s creates a Bloch electron with combined band and spin index s at wave- vector k, 
and an electron-electron interaction term: 

U = \ ^q,s,s'4+q, S C k'-q,s'Ck', S 'C k , S , (A2) 

k,s,k',s',q^0 

where Uq, s ,s' is the Coulomb matrix element. 

The exciton operator B n u of a given transition a is associated with the relative motion 
wave function <p n ^ a at zero center-of-mass momentum: 

B n ,a = Ck, CT c kjS ' c k, s - (A3) 

k 

The transition a = a(s, s') connects an electron from a valence band with combined band 
and spin label s' to a conduction band state with label s. The corresponding ^azr-operator 
c iks' c k,s is denoted in Sec. |TJ as Vk,<r- Optical selection rules determine whether the a- 
transition is an optically allowed dipole transition with matrix element /x CT or a so-called 
dark transition. Dark states are connected to optically active states via a spin-flip process, 
which is assumed to be very slow on the time scales of interest. The selection rules also 
determine the corresponding helicity of the dipole transition. 

The first step leading to the force operator D is the commutator C n ^ = [B n>a , H], 

Cn,a = XX e M 2 - ek, Sl )0*, ki(T c kiS2 c kiSl (A4) 

k 

~~ ^q,si,S2 0n,k,<T C k+q,S2 C k+q,si + ^q,si,si0n,k,cr C k,S2 C k,si 

k.q^O k,q^0 

+ Y (^q,s,s 2 0n,k'-q,o- ~~ Uq,s,si 0n,k',cr) C k+q,s C k,sC k /_ qSl Ck', S2 • 
k,k',s,q^0 

Using c ks Ck,s2 = En^fik the Wannier equation for the exciton wave function simpli- 

fies the first three terms on the rhs of Eq. (A3): 
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k,k' 



ek,s 2 — e k,si + X! ^q,si>si ) ^k,k' — ^k-k',si,s 2 ( 1 ~~ ^k,k' 
q^O 



^k.n.o-^k'/ra,^ 



(A5) 



and, therefore, 



B n>a (A6) 

+ ^2 (^q,s,S2 0n,k'-q,o- — ^1,s,S! 0n,k',o-) C k+q,s C k,s c k'-q,si C k',s 2 • 
k,k',s,q^0 

The force operator is given by the commutator D p j. n ^ = [B p j, C n>er ]. The first contribu- 
tion on the rhs of Eq. ( |A6[ ) vanishes because [B Pt i, B n>a \ = 0. Thus, 

D p j in ,a = (^q^2,s 2 0n,k'-q,cr ~~ ^q,i2,si 0n,k',<r) ^p,k+q,Z 

k,k',q^0 

— (^q,/l,s 2 < / , n,k'-q,<7 ~~ ^q,«l,si0n,k',<r) ^p,k>* 



X C k;+q,«i C k,«2 C k'-q,si C k',S2 



(A7) 



The spin-independent Coulomb interaction leads to a further simplification. In terms of the 
operator 



Ai,<r(q) — X! (0n,k-q/2,<r ~~ 0n,k+q/2,<r) C k-q/2,si C k+q/2, 



',«2 > 



(A8) 



the .D-operator can be written in a compact form: 

Dp,l;n,a = ^ UqApj (q)A»,<T (~q) • 
q^O 



(A9) 



The A-operator can be related to excitons with finite center-of-mass momentum Q. The 
mass ratio of electrons and holes play a crucial in the correlation function dynamics. We 
define for the transition a = (s, s') the positive mass-ratio K a = m s /m s r. With 



V* 



1 + K n 



(A10) 



the generalization of Eq. (|A3|) reads 



B n ,a(Q) — ^n ) k,cr c k-i7 CT Q,s /C k+K £ r7? CT Q 



(All) 



We have neglected a possible Q-dependence of the relative-motion wave-function, which is 
valid for parabolic bands. Transforming the eh representation in Eq. flA~8f) to one in terms 
of the excitons leads to 



^n,a(q) = (</>n,k-r, CT q,<x _ C,k+7, CT K CT q, CT ) <j>a,V,aB a>a {<\) 
k,Q 

= E/n CT a(q)5a, CT (q). 



(A12) 
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The -D-operator in Eq. 



becomes in terms of the excitons: 



«,a',q^0 



(A13) 



We use £/ q = t/- q , to ensure the symmetry of the D-operator in the exciton labels. 

The initial value r = of the correlation function Eq. Q4.1ij| ) is a ground-state (vacuum) 
expectation value. We first calculate 



+ (0|S aV (q) \[B a!a (-q),Bl )a „(-q% B^W] l°>- 



(A14) 



The initial value 7^^'" ' CT of Eq. (|4.13 ) can be split in three contributions. The first term 
on the rhs of Eq. (|A14 ) gives 



7 « = 5^,5^,, y: ^/!^(q)/nV(-q)/n^(q)7n^(-q)*- 

A/3',q 



(A15) 



The summation over the exciton labels /3, f3' can be performed, using the abbreviation 



<,Jq) = E/nVq)A/3(q) 



(A16) 



which gives 



7 (1) = ^, CT '^, CT "E^q^,n"(q)<n'(~q) 



(A17) 



and for the second term on the rhs of Eq. ( |A14 ) 

7 {2) = V'^EteWC"! - ^ 



(Ai8) 



The third contribution arises only if the two electron-hole pairs have at least one band in 
common, i.e., in the parallel-spin or coupled-spin case, c.f. Sec. ITJ. In the case of excitons 
being ideal bosons, only 7W and 7^ would contribute. We give the result for the general 
case: 



-,( 3 ) — ( a( c ) a( c ) a^) a^ i a'"' a^) a( c ) a( c ) \ tt tt 

T - ~ {°a,a" d v,a> d a,a> d a,a" + 6 o,a" 6 a,a' d <J,o' 6 a,o" ) 1^ ^q^k-k' 



(A19) 



q,k,k' 



X 



y$ft,k'-q/2,CT ~ $ ft,k'+q/2,s)(*n",k'-q/2,<7" ~ $n",k-q/2,<X 



X ( < &n',k+q/2, CT ' - ( &n',k'+q/2, ( 7')( $ n,k+q/2, ( T _ $ n,k-q/2 l( r) • 

The Kronecker-delta symbol with an upper index means, the corresponding conduction (c) 
or valence-band (v) must be identical for the transition pair. 
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APPENDIX B: AN INTERACTING SEMICONDUCTOR MODEL 



We present the details of the real-space extended Hubbard-model used for the numerical 
calculations of the correlation-function and the linear optical properties (excitons). We start 
with the general <i-dimensional lattice-model of an electron-hole system and briefly review 
the fundamental aspects of the problem on a lattice. The kinetic energy part of a two-band 
lattice model is the usual hopping term 

#0= E t^ m/ cl tS c mf , s + h.c, (Bl) 

m,m',s 

where s creates an electron in a Wannier-state at the lattice site m with a combined band 
and spin index s and ^' m , is the hopping matrix element between two sites. The hopping is 
restricted to nearest neighbor with parameter for each band. We take < (£W > 0) 
for an s-type (p-type) conduction (valence) band, allowing for interband optical transitions. 
The potential energy of the system is given by the Coulomb interaction between the electrons 
and charge-positive lattice ions with charge eZ core > 0. In terms of the dimensionless charge- 
density (electrons and ions) at site m 



Pm ^^-^ C m,s C ni,s^ ^core > (B2) 

the Coulomb interaction in the Wannier-states to leading order is given by the electrostatic 
monopole-monopole contribution 

U = - E ^m.m'PmPm'- (B3) 
m,m' 

where U m , m i = e 2 /|m— m'| for different sites and the on-site Coulomb interaction £/ mjm = Uq 
is an additional parameter in the model. 

The ground-state |0) of the lattice model is given by the completely filled valence-band 
states at each site. For the numerical calculation of the correlation function, a four-band 
approximation of the semiconductor is sufficient, i.e., only two holes and two electrons of 
one or two species are relevant. The index s is split into two bands and two spin directions. 
This leads to different types of correlation functions, as discussed in Sec. |II|. We write the 
ground state as 

\0) = U m cl^cl A Jvac). (B4) 

The spatial extension of the relative motion of a single exciton in its ground state depends on 
the ratio of the sum of the bandwidths B s = B\ + B 2 = 2d(i 1 — 1 2 ) and the nearest neighbor 
Coulomb energy U\ = e 2 /a L , where is the lattice constant. For B s /U\ ^> 1 the results of 
the lattice model are similar to the continuum limit given by the usual two-parabolic-band 
effective-mass model. If one expands the band dispersion quadratically around the center of 
the Brillouin zone one finds for the Bohr radius of the exciton for d = 3 

2(ti-t 2 ) 

a = a L (B5) 
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and the exciton binding energy lo x is given by 



■* = £(^V (B6) 



2 \a J 

In the following, we give all energies in units of the Rydberg energy u x . One should keep 
in mind that this is not the exciton binding energy in d = 1, which apart from a different 
dimensionless prefactor depends on the ratio rj = Uq/U\. The different effective electron 
(hole) masses enter via |ti/£2| — I^A^il- The total Hamiltonian operator H = Ho + U 
is conveniently transformed into the electron-hole picture defining c m)S = c m ^,s for electron 
states and h m ^ s = cl n l s for the hole states. Let |$) denote a many-particle eigenstate of H 
with an arbitrary number N p of the eh-pairs as a quantum-number. As the corresponding 
pair-number operator N p commutes with H, we can always work with eigenstates of N p . This 
is the fundamental assumption for the Hubbard-operator formulation of the Hamiltonian in 
Sec. 0. 

We now discuss in detail the results of the numerical calculations for a quasi-one- 
dimensional ring model. The Coulomb interaction Eq. (|B3|) between two charges depends 



on the chord-distance 161 



|m — m | = sm[— |m — m \\ (B7) 



7T 



of the sites, which are labeled by the dimensionless numbers m. The model interpolates 
smoothly between the two limits of Frenkel excitons cto/ a L ~~ > and the Wannier limit of 
large extended objects cto/ a L ^ 1, which is only limited by the finite total system size N. 
The system size for the numerical calculation of the correlation function is N = 120 for 
equal-spin correlation and iV = 140 for opposite-spin correlation. The on-site Coulomb 
interaction is fixed with 7] = 1.5. We have also fixed the bandwidth B := |4(t2 — t\)\ of 
the eh-pair continuum for each model, which depends only on the reduced electron and 
hole mass, while the positive mass ratio m e /mh is varied. The mass ratio is an important 
parameter of the two eh-pair subspace. We compare the positronium limit m e /mh = 1 
with the semiconductor GaAs case of m e /mh ~ 0.15 and the molecular (hydrogen) limit 
m e /mh 1. The long-range Coulomb potential leads to the formation of a finite number 
of bound exciton-states in the system in contrast to the usual Hubbard-model with purely 
on-site repulsion U . 



1. Parallel-spin excitons 

The parallel-spin exciton case is characterized by the exchange repulsion between same 
type carriers in the bands, which, in an ideal one-dimensional system, plays the role of dy- 
namical boundary conditions. The parallel-spin case is relevant when the optical excitation 
process is limited to a single circularly polarization of the external field. A complete set of 
parallel-spin two-pair states with zero center of mass momentum is given by pTl 



\p,a,f3) = y= ^ cL+p+afyL+p C TO.+fl^m 1 0) J (B8) 

V i V 777, 
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where 

(j(a,P) _ f p = N/2 and a = (3 , Bg , 

p | 1 e/se 

is a normalization constant and a triplet (p, a, (3) labels the relative position of the carriers. 
In order to avoid double counting of states, the set of possible triplets is restricted. The 
(single) spin-index is understood. 

The real-space representation of the product state of two excitons with quantum numbers 
a and b with band indices (a 2 , ai) and (b 2 , which enters the calculation of the generating 
correlation function, is given by 

M lb\°) = Tf E fe (m 1 -m 2 )0 a (m 3 -m 4 ) (BIO) 

mi ,7-712, 1713, rrt4 
X C mi,6i^m2,fe2 C m3,ai^ , m.4.a2 1^)' 

where each exciton has zero center-of-mass momentum. The real-space exciton wave- 
function is given by 

\ k = -^'E e ~ ikm M™), (bii) 



and the exciton operator can be expressed as 

B\ = <Pa{mi-m 2 )c ] miAl h ] m2A2 . (B12) 

V-'' m\,rri2 

For the force-force correlation function in the nonlinear optical response, we need a linear 
superposition of two-pair states in the initial state Dl b \0) 

D l,b\°) = T? E <M m i ~ m 2 )(p a (m 3 - m 4 ) (B13) 

7711,7(7,2,7713,7774 

X (Urns— mi U m4 — mi U m ^— m2 ~\~ U m ^—m-z) 
X C mi,6i^m2,()2 C m3,ai ^7114,0,2 1^) ' 

Both states are center-of-mass eigenstates with zero total momentum in systems with peri- 
odic boundary conditions. It is possible to fix the hole position m 4 and to introduce relative 
distances with respect to this hole to reduce the number of coordinates. The states Eq. ( |B10| ) 
and Eq. ( B13|) can then be expressed in the basis set Eq. (|B~8[). 



2. Opposite-spin excitons 

The corresponding basis set for the opposite-spin pair-states is chosen as 

|mi, 777,2,777,3) = —j= E C m 3 +m4,a 1 hln2+m4,a2 C ln 1 +m4M^m4,b 2 \^)^ (B14) 

ViV 7714 
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where the opposite-spin indices for the carriers in one band are understood, i.e., 62 7^ 0,2 and 
bi 7^ a±. For numerical work, the absence of the Pauli-blocking for the two electrons/ holes 
in the bands leads to a quite simple counting of opposite-spin states in contrast to the 
equal-spin case. 

The numerical evaluation proceeds in the usual way. The total number of basis states 
involved in the calculation is N 3 for the opposite-spin problem, which gives a vector length 
of about 3 x 10 6 . The Lanczos algorithm requires two states for the iteration, which reside 
in memory to speed up the process and can be handled quite well on a PC with a total of 64 
MByte memory. The Lanczos algorithm tridiagonalizes the Hamiltonian for the 4-particle 
problem, starting with the (normalized) initial state D^O) for the force-force correlation 
function. Each iteration step produces one new basis state. The iteration is extremely 
fast, due to the sparseness of the Hamiltonian matrix of H in real-space. We truncate the 
iteration after the spectrum of the resolvent matrix (z — H)^ 1 stabilizes. No eigenvectors or 
eigenvalues for the biexciton problem have to be calculated. The relevant spectrum is given 
in the usual way 

U(uj) = -2Im(0\D a j b u _* +iQ Dlj b \0) (B15) 

from the inversion of the resolvent matrix, which is simple in the tridiagonal representation 
of if. 
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FIGURES 



FIG. 1. Selection rules of heavy and light-hole exciton transitions in Zincblende structures. 
A splitting of heavy- and light-hole states in heterostructures allows for a selected excitation of 
these transitions (I). The degenerate case (II) is shown for each helicity of the exciting field. The 
A-transition (III) is excited by linear polarization. 

FIG. 2. Force-force correlation function spectra F(u) for the ls-exciton contribution for op- 
posite spins (solid line) and parallel spins (dashed line). The mass-ratio of electrons and holes 
is m e /mh = 1 for (a), which corresponds to the positronium limit, m e /mh = 0.15 for (b), which 
corresponds to the heavy-hole/electron mass ratio of GaAs and for the molecular (hydrogen) limit 
i r ne/ m h = in (c). Bound excitonic molecules appear for u < and continuum two-exciton 
contribution have uj > 0. 

FIG. 3. Three-pulse four- wave- mixing transmission geometry. 

FIG. 4. Source term C(t, 0) in the nonlinear response in the parallel-spin case with parameters 
of Fig. ||). The mean- field approximation, which has no real part from Eq. ( |4.14| ) overestimates 
the response. 

FIG. 5. Source term C(t, 0) for the nonlinear response in the opposite-spin case with the spec- 
tral function of Fig. |2| (b). The dephasing time is T2 = 2 ps. The pronounced oscillations with a 
period of 2.8 ps {E xx = 1.5 meV) origin from the bound-state biexciton in the system. 

FIG. 6. TR-signal of the nonlinear polarization for the one-dimensional semiconductor model 
for T = in the weak dephasing limit. The oscillations with the biexciton binding energy are 
a ringing in the signal, since no additional biexciton states are necessary for the response. The 
dephasing time is T2 = 4 ps. 

FIG. 7. TR-signal of the nonlinear polarization for the one-dimensional semiconductor model 
for T = in the strong dephasing limit in comparison with the mean-field response for parallel-spin 
excitation. The dephasing time is T<i = 0.5 ps. 

FIG. 8. Tl-intensity of the FWM-signal for co-polarized (solid line) and cross-polarized (dashed 
line) circular excitation for a dephasing time of T2 = 4 ps. 

FIG. 9. Normalized Tl-intensity of the FWM-signal for co-polarized circular excitation on a 
log scale. For negative time delay, significant deviation from the exponential decay with a rise-time 
of T2/4 is observed for smaller dephasing. 
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FIG. 10. Normalized Tl-intensity of the FWM-signal for cross-polarized circular excitation on 
a log scale. For negative delay times, oscillations with the binding energy of the biexciton are 
visible for sufficiently small dephasing. 
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